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Abstract 

We present results of a lattice analysis of the B parameter, Bb, the decay constant 
/b, and several mass splittings using the static approximation. Results were obtained 
for 60 quenched gauge configurations computed at (3 = 6.2 on a lattice size of 24 3 x 48. 
Light quark propagators were calculated using the 0(a)-improved Sheikholeslami- 
Wohlert action. We find B B taXlc (nib) = 0.69 + Jj(stat) + ^(syst), corresponding to 

B| tatic = L02 +5 +3, and /|tatic = 26 g +W +* M eV, fa B B J ft B B = 1M ^ 

where a variational fitting technique was used to extract Jj^ atlc . For the mass split- 
tings we obtain M Ba - M Bd = 87 ^MeV, M Ab - M Bd = 420 + 1 ™ ^MeV and 

Ml*-M% = 0.281 ^GeV 2 . We compare different smearing techniques intended 
to improve the signal/noise ratio. Prom a detailed assessment of systematic effects we 
conclude that the main systematic uncertainties are associated with the renormalisa- 
tion constants relating a lattice matrix element to its continuum counterpart. The 
dependence of our findings on lattice artefacts is to be investigated in the future. 
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1 Introduction 



Heavy quark systems have attracted considerable interest in recent years. Studying the de- 
cays of hadrons containing heavy quarks is important in order to narrow the constraints on 
the less known elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Precise knowl- 
edge of the CKM matrix elements serves to test the consistency of the Standard Model and 
to detect possible signals of "new physics". Theoretical tools for dealing with heavy quark 
systems, such as the Heavy Quark Effective Theory (HQET) [|IJ 0, [|, have been developed 
and are being successfully applied in the analysis and interpretation of experimental data. 
However, theoretical estimates of form factors, decay constants and mixing parameters are 
subject to uncertainties due to strong interaction effects whose nature is intrinsically non- 
perturbative on typical hadronic scales. Lattice simulations of QCD are designed to provide 
a non-perturbative treatment of hadronic processes and have already made important con- 
tributions to the study of the spectroscopy and decays of hadrons in both the light and heavy 
quark sector Q. For systems involving heavy quarks, most notably the b quark, the role 
of lattice simulations is two-fold: firstly, to make predictions for yet unmeasured quantities 
such as the decay constant of the B meson, fs, or the masses of baryons containing b quarks; 
secondly, to test the validity of other theoretical methods such as large mass expansions or 
the HQET. 

One problem that is encountered in current simulations of heavy quark systems is the fact 
that typical values of the inverse lattice spacing lie in the range 2 — 3.5 GeV which is well 
below the b mass. There are several methods for dealing with this problem, one of which was 
proposed by Eichten || in which the heavy quark propagator is expanded in inverse powers 
of the heavy quark mass. The so-called static approximation is the leading term in this 
expansion, for which the b quark is infinitely heavy. One may also hope to compute some of 
the higher-order corrections to the static limit, although the presence of power divergences 
presents theoretical and practical complications 0. 

Another method for lattice studies of heavy quark systems is to use propagating heavy 
quarks. At present, these simulations are carried out for quark masses around the charm 
quark mass, and the results obtained in this fashion must be extrapolated to the mass of 
the b quark. Clearly, this method depends crucially on controlling the effects of non-zero 
lattice spacing ("lattice artefacts") at the heavy masses used in the simulation. In general, 
the influence of lattice artefacts on quantities involving propagating quarks can be reduced 
by considering "improved" actions as suggested by Symanzik @ and detailed further by the 
authors of §| and |J. For heavy-light decay constants, improvement has been successfully 
applied to quark masses in the region of that of the charm quark [|TU], |TTJ. Furthermore, 
the data from the static approximation, obtained at infinite quark mass, serve to guide the 
extrapolation of results obtained using propagating heavy quarks to the mass of the b. 



In this paper we report on our results for f B , the B parameter B B describing B° — B° 
mixing and mass splittings involving the B and B* mesons as well as the The results are 
obtained using the static approximation for the heavy quark, whereas the 0(a)-improved 
Sheikholeslami-Wohlert action || [| is used for the light quarks. In many ways this study 
is intended as a continuation and extension of earlier simulations. For example, we are able 
to study the effects of 0(a)-improvement on various quantities, among them the B* — B 
mass splitting, which is believed to be sensitive to discretisation errors. Furthermore, this 
splitting arises only at next-to-leading order in the large mass expansion and serves therefore 
as a measure of higher-order corrections to the static limit. 

We now list our main results. For the B parameter at scale m& in the static approximation 
we obtain 

B*f c (m b ) = 0.69 + _l (stat) + 2 (syst), (1) 
which corresponds to a value of the renormalisation-group-invariant B parameter of 

Bgf = 1.02 ! 5 6 (stat) ^(syst). (2) 

The decay constant Jj* atlc is found to be 

^tatic = 26g +18 (stat) +28 (sygt) MeV> (3) 

From the chiral behaviour of _gg atlc and /^ tatlc we extract 

=1.34 +;(stat) ^(syst). (4) 

J B d B d 

Finally, for the mass splittings we obtain 

= 87 (stat) +^(syst)MeV (5) 
= 420 + 1 ™ (stat) ^ (syst) MeV (6) 
= 0.281 (stat) +f 7 (syst) GeV 2 . (7) 

Here, the systematic error on dimensionful quantities is dominated by uncertainties in the 
lattice scale, whereas systematic errors on dimensionless quantities arise from the spread of 
values using alternative fitting procedures. There is an additional uncertainty in the estimate 
of the B parameter, which arises from the perturbative matching between full QCD and the 
lattice theory in the static approximation. We estimate this uncertainty to be as large as 
15-20%, as will be discussed later. 

Our estimate for the B* — B hyperfme splitting in eq. (0) is much smaller than the experimen- 
tal value of 0.488(6) GeV 2 (as is also the case with Wilson fermions The implications 
of this result are discussed in more detail in section |4~3| . Here we only wish to state that the 
relevant matching factor for the B* — B splitting is also subject to significant uncertainties. 
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It is beyond the scope of this paper to investigate the dependence of our results on the 
lattice spacing a. Here we wish to stress that we seek to reduce these effects considerably 
by employing an 0(a)-improved fermionic action, leaving an extrapolation to the continuum 
limit to future studies. 

The paper is organised as follows. In section |2] we describe the details of our simulation 
and analysis. Section |3] contains our results for _gg atlc and J^ tatlc . The splittings B s — B^, 
B* — B and — B are discussed in section [|. Finally, section |5| contains a summary and 
our conclusions. 



2 The Simulation 

Our results are based on 60 SU(3) gauge configurations in the quenched approximation, 
calculated on a lattice of size 24 3 x 48 at f3 = 6.2. The configurations were generated using 
the hybrid over-relaxed algorithm described in O . 



Light quark propagators were computed at three values of the hopping parameter Ki, namely 
k, = 0.14144, 0.14226 and 0.14262, using the Sheikholeslami-Wohlert (SW) action § 

S S F W = Sj-i%Y, F^x)^ q(x), (8) 

X,fJ,,U 

where Sp is the standard Wilson action and F^ v is a lattice definition of the field strength 



tensor |TJj 



The chosen Ki values are in the region of the strange quark, whose mass, as determined from 
the mass ratio m 2 K /m 2 p , corresponds to k s = 0.1419(1), while the chiral limit of massless 
quarks is reached at n crit = 0.14315(2) fl~5|j . 

The leading term in the expansion of the heavy quark propagator is given by 

S Q (x,t;0,0) = |e(t)e- m ^^^ + e(-t)e m «'i^| 8®(x)V s (t,0), (9) 

where V$(t, 0) is the product of links from (0, t) to the origin, for example for t > 0, 

V s (t, 0) = Ut(0, t - l)Ut(0, t - 2) • - • [4(0, 0). (10) 

The static quark propagator, eq. ([5J), and the light quark propagators were combined to 
compute correlation functions for the relevant quantities in this paper. 

In order to obtain 0(a)-improved matrix elements for heavy-light bilinears using the static 
approximation, it is sufficient to carry out the improvement prescription in the light quark 



sector only [27]. Therefore we describe the light quark using the SW action and consider 
operators in which only the light quark field q(x) has been "rotated" H |10j, i-e. 



O v = b\x)T(l-\ T D)q(x). (11) 



Here, b(x) denotes the heavy quark spinor in the static approximation, and T is some Dirac 
matrix. 

It is useful to use extended (or "smeared") interpolating operators in order to isolate the 
ground state in correlation functions efficiently. This is of particular importance in the static 



approximation where calculations using purely local operators are known to fail JT6[ . In this 
study we compare different smearing techniques which can be broadly divided into gauge- 
invariant and gauge-dependent methods, the latter being performed on gauge configurations 
transformed to the Coulomb gauge. 



For gauge-invariant smearing we use the Jacobi smearing algorithm described in [jiT]] . The 
smeared field at timeslice t, q (x, t) is defined by 

g 5 (f,£) J(x,x') <?(£', t) (12) 

x' 

where 

N 

J{x,x') = J2 «sA"(x,a/), (13) 

n=0 

and A is the gauge-invariant discretised version of the three-dimensional Laplace operator. 
The parameter k$ and the number of iterations iV can be used to control the smearing 
radius. Here, we quote our results for k$ = 0.25 and N = 140, which corresponds to a r.m.s. 



smearing radius of r = 6.4 |T(J. The same values were used in our previous study of f§ 



static 

j 



obtained on a subset of 20 of our 60 gauge configurations [T(J 



In order to study smearing methods in a fixed gauge, the configurations were transformed 



to the Coulomb gauge using the algorithms described in |L8], |19| . The lattice analog of the 
continuum Coulomb gauge condition, diAi{x) = 0, is 

9{x) = Tr(D(x)D\x)) = (14) 

where 

D(») = i;(Oi(x) + D?(x-i) - h.c.) t _ ta (15) 

i 

with the index i signifying spatial components only. At each iteration of the gauge fixing 
procedure the average value of 9 was calculated, (9) = J2x@( x )/V, where V is the total 
lattice volume. For each gauge configuration the gauge was fixed to a precision (9) ~ 10 -4 . 

Defining the smeared quark field q s (x, t) via 

q s (x,t) = ^2f^^')q(x',t) (16) 

g> 

we considered the following smearing functions f(x,x') for smearing radius r : 



Exponential: f(x,x') = exp { — \x — x'\/r } 



(17) 



Gaussian: f(x,5?) = exp |— \x — a 7 '] 2 / 7 "^} 



3 



Cube : f(x, x') = J\ 6(r - \x, t - x[\) (19) 



i=l 

Double Cube: f(2,£) = f[ fl - ^ ~ ^ ) 6(2r - k - x'A). (20) 

Following the analysis in ref. []2U|] , where a variety of smearing radii was studied, we chose 
r = 5 in all cases. 

Our 60 gauge configurations and the light quark propagators were computed on the 64- 
node Meiko i860 Computing Surface at Edinburgh. The transformation to the Coulomb 
gauge was performed on the Cray Y-MP/8 at the Daresbury Rutherford Appleton Labora- 
tory. Smeared propagators using the gauge-invariant Jacobi algorithm were calculated on a 
Thinking Machines CM-200 at the University of Edinburgh. All other smearing types and 
the relevant correlators were computed on a variety of DEC ALPHA machines. 



Statistical errors are estimated from a bootstrap procedure |2T) , which involves the creation 
of 200 bootstrap samples from our set of 60 configurations. Correlators are fitted for each 
sample by minimising x 2 , taking correlations among different timeslices into account. The 
quoted statistical errors are obtained from the central 68 % of the corresponding bootstrap 



distribution |L4| 



In order to convert our values for decay constants and mass splittings into physical units we 
need an estimate of the inverse lattice spacing in GeV. In this study we take 



-i 



2.9 ±0.2 GeV, (21) 



thus deviating slightly from some of our earlier papers where we quoted 2.7 GeV as the central 
value [ 1(| [14], . The error in eq. ([H]) is large enough to encompass all our estimates for 
a -1 from quantities such as m p , f w , tun, the string tension \f~K and the hadronic scale Rq 



discussed in |2"2]] . The shift was partly motivated by a recent study using newly generated 



UKQCD data p3[: using the quantity R we found a 1 = 2.95 ^ GeV. Also, a non- 



perturbative determination of the renormalisation constant of the axial current resulted in 
a value of Za = 1.04 [[24]] which is larger by about 8 % than the perturbative value which we 
had used previously. Thus the scale as estimated from f n decreases to around 3.1 GeV which 
enables us to significantly reduce the upper uncertainty on our final value of a -1 [GeV]. 



3 Decay Constants and Mixing Parameters 

In this section we present the results for Jj^ tatlc and I?§ atlc . We begin by listing the various 
operators in the lattice effective theory and discussing the relevant renormalisation factors. 



Here, we wish to emphasise that all our matching and scaling factors are consistently defined 
at leading order in the strong coupling constant. The 2- and 3-point correlators are defined 
before the results are discussed. We close the section with a discussion of the phenomeno- 
logical implications of our findings. 



3.1 Lattice Operators and Renormalisation 

In the continuum full theory, the pseudoscalar decay constant of the B meson is defined 
through the matrix element of the axial current: 

(0\A,(0)\B(p))=tp,f B , A^{x)=b{x) llllb q{x). (22) 

The AB = 2 four-fermi operator Ol which gives rise to B° — B° mixing is given by 

L = (67„(1 - 75) q) - 75) q). (23) 

This operator enters the effective Hamiltonian describing B° — B° mixing whose amplitude 
is usually expressed in terms of the B parameter, which is the ratio of the operator matrix 
element to its value in the vacuum insertion approximation 



o\ 



b bw - — 8 72 M 2 — ' y 2A ) 

3 JB m B 

where /i is a renormalisation scale. We have adopted a convention in which f n = 132 MeV. 
The scale dependence of B B (/J>) can be removed by multiplication with a factor derived 
from the anomalous dimension of the operator Ol- At one-loop order, one can a define a 
renormalisation-group-invariant B parameter by 

B B = a s (fxy 2 ^°B B (fi), (25) 

where fi = 11 — 2n//3. The strong coupling constant appearing in the above expression is 
usually defined in the MS scheme. We will use the one-loop expression for B B below, since it 
is consistent with our matching between lattice and continuum results, which is performed 
at order a s . Alternative matching procedures and higher-order scalings of B B (fj,) can always 
be incorporated by suitably replacing our matching and scaling factors. 

In order to get estimates for the matrix elements of the axial current and of the four-fermi 
operator Ol in the continuum, these operators need to be matched to the relevant lattice 
operators in the static effective theory. The matching of operators in the static approximation 
is usually performed as a two-step process, in which one first matches the operators in the 
effective theory on the lattice to those in the continuum effective theory. In the second step, 
one then matches the continuum effective theory to the operator in full QCD. 



For the axial current, this two-step matching process was considered for Wilson fermions to 
one-loop order in [ 25| , |26| and extended to the 0(a)-improved case in ]27[ and |28| . 



At n = a 1 the renormalisation factor between full QCD and the lattice effective theory at 
order a s for the SW action is 



alia x ) r , o ai att fa _:l 



^static = 1 + ^_J_ | 2 ln(G 2 m 2 } _ 2 j _ 15 02 " s J (26) 

In order to evaluate Z^ atlc numerically for m& = 5 GeV, we define the strong coupling 
constant for rif active quark flavours at leading order in the continuum by 

A) hi (/i/A^O 

where & = 11 - |n/, and we take n f = 4, = 200 MeV, fj, = a' 1 = 2.9 GeV. Thus, 
despite the fact that our results for matrix elements of lattice operators are obtained in 
the quenched approximation, we use the relevant number of active quark flavours at the 
respective renormalisation scale when matching the continuum effective theory to full QCD. 
This concept implies that we abandon the quenched approximation once we have obtained 
the matrix elements in the continuum effective theory after the first matching step. Of course, 
all our results are still subject to a systematic error due to quenching, which is, however, 
hard to quantify unless precision data from dynamical simulations become available. 

For the matching step between the effective theories in the continuum and on the lattice we 



take the "boosted" value for the gauge coupling [29, 30 



where u$ is a measure of the average link variable, for which we take uq = (S^cnt) -1 with 
K crit = 0.14315(2) |]|. With these definitions, we find 

^static = q 7 g ( 2 g) 

This is very close to the value of Z^ atlc = 0.79 quoted in our previous study |[10|| , and also to 
^static _ o gi use d in a simulation by the APE collaboration [36| employing the SW action 
for light quarks at (3 = 6.2. 



In the case of the four-fermi operator the situation is more complicated due to operator 
mixing. When relating the full continuum theory to the continuum effective theory, it is 
useful to introduce 

O s =(b(l- l5 )q)(b(l~ l5 )q). (30) 

This operator is generated at order a s in the continuum owing to the mass of the heavy 
quark |3T|. The one-loop matching factors between the continuum full theory at scale m& 



and the continuum effective theory at scale /i < m;, are given by 4 



OT\m h ) 



with C L = -14 and C & 



167T 2 
3 



4 In K//. 2 ) + C L ] Of 0*) + ^O* Of (?) (31) 



In the matching step between the continuum effective and the lattice effective theories, two 
additional operators are generated due to the explicit chiral symmetry breaking induced by 
the Wilson term 



Or = (b lll {l + l5 )q)(b lll {l+ l5 )q 

N = (67/i(l-75)9) (&7/*(l + 75)g 

+ (b lll (l + l5 )q)(b llI (l- l5 )q) 

+ 2(fo(l-7 5 )g)(6(l + 75 )g) 

+ 2(6(l + T8 )<z)(6(l-7 B )g). 



(32) 



(33) 

For the 0(a)-improved SW action, the full one- loop matching for the four-fermi operator 
Of to the lattice effective theory at scale \i = a -1 is given by 



O 



cS/ 




47T 



[D R + D R ] 0>f + 



4:71 



[D 



N 



D N ] 0\ att 



'N 



(34) 



The coefficients D Lj D R , D N were calculated in |3l[ , whereas those for the SW action, 
D ! Ll D R , D ! N are listed in [p7| . A subtle point to note is that the coefficient D L is quoted 
as D L - 



-65.5 in refs. 0, 0. In it was stated, however, that the reduced value of 
the quark self-energy should be used if the static-light meson propagator is being fitted to 
the usual exponential. Using the reduced value in the formula for yields a value of 
Dl = —38.9, and hence results in a much smaller correction to Of in the matching step 
between the lattice effective and the continuum effective theory. In the following we quote 
numerical values for all relevant renormalisation constants using the reduced value of the 



quark self-energy. It should be added that the expression for i n ec L- (PBj ) is also based 

on the reduced value , and thus our procedure to extract the B parameter from a suitable 
ratio of matrix elements in the static theory is entirely consistent. 

Expanding the various matching factors to order a s , we arrive at the following expression 
for the matching of the operator 0^ u (m&) to the operators in the lattice effective theory: 



oi uii ( 



1 + 



aAa 



An 



41n(a 2 m 2 ) - 14 



- 22.06 



47T 



O 



latt 



Note that in refs. |3lL M the operator is obtained at the scale \x 



2 < mi,. This requires 

the factor 4 multiplying ln(rn^//^ 2 ) to be replaced by a factor 6, which is the difference of the anomalous 
dimensions in the continuum full and continuum effective theories. 



_ 4.19 ^ O 1 ^ - 13.96 ^ Off* - 2 Q^t (35) 

= ZlO 1 ^ + Z R 1 ^ + Z N X ^ + Z s Of\ (36) 

It is this expression which we will use from now on to convert our lattice results into an 
estimate of the B parameter. 

Using our numerical values for the coupling constants a s (a~ l ) and a^ att (a _1 ) we find 

Z L = 0.55, Z K =-0.04, ^ = -0.15, Z s = -0.18. (37) 

As was already mentioned in |TT|, the correction to the matrix element of 0^ tt is rather large, 
thus calling the applicability of one- loop perturbative matching into question. In fact, if the 
matching is performed by first computing Of^a -1 ) according to eq. (|34|), and then inserting 
the result into eq. (|31|), we observe that our estimate for Bb increases by 20%. This way of 
matching includes some part of the 0{a 2 s ) contributions to the renormalisation factors, and 
hence it serves to estimate the influence of higher loop corrections in the matching procedure. 



In reference [g2j the two-loop anomalous dimensions of the axial current and the four-fermi 
operator were calculated for the effective theory in the continuum. Including this result into 
the matching step between 0^ u (mb) and Of* (//) in eq. (|3lD changes the final result only by 
1-2 %. Therefore we conclude that the bulk of the uncertainty arises from the matching step 
between the continuum effective and lattice effective theories, and also from the large factor 



of Cl = — 14 in eq. (31). An entirely non-perturbative determination of the renormalisation 
constants relating Off to the different lattice operators, as outlined in [33|, is highly desirable 
in order to clarify this important issue, which is of equal importance in the case of Z^ atlc , as 
will be illustrated later. 



3.2 Correlators for 2- and 3-point Functions 

In order to extract the pseudoscalar decay constant we consider the correlation function of 
the time-component of the improved static-light axial current 

£<^(x,t)4(0,0)> ^ &f^e- M *\ (38) 

X 

In practice, using particular combinations of the smeared (S) and local (L) axial current, we 
compute correlation functions defined by 

C ss (t) = ^(0\A s 4 (x,t)A\ S (6,0)\0) *a (Z s ) 2 e- Et (39) 

X 

C LS (t) = ]T(0|^(£,t)4 5 (0,0)|0> ^Z s Z L e- Et , (40) 

X 

C SL (t) = ]T(0|Af(£,t)4 L (0,0)|0> ^Z¥e- E1 , (41) 

X 



where E is the unphysical difference between the mass of the meson and the bare mass of 
the heavy quark. 

The pseudoscalar decay constant yj* atlc is related to Z L via 



f B tatic = Z L ^jj-ZT ic , (42) 
where Mb is the mass of the B meson. 

Z L and thus /^ tatlc is extracted from C ss (t) and C LS (t) as follows: by fitting C ss (t) to the 
functional form given in eq. fl39|) we obtain Z s and E. At sufficiently large times the ratio 
C LS {i)/C ss {t) — > Z L jZ s , so that Z L can be determined. As was observed earlier [j34[, |35|, [TT 



using the correlation function C LS (t) in which the operator at the source is smeared yields 
much better statistics than the corresponding correlator C SL (t) for which the smearing is 
performed at the sink. 

Alternative methods, discussed in [SB], include a direct fit of C LS {t) in order to extract Z L Z S 



which can then be combined with the ratio Z L jZ s to compute Z L . However, this method 
requires the ground state to be unambiguously isolated, which is more difficult for C LS (t), 
since the plateau in the effective mass plot is approached from below. 

A more direct method, which does not involve any fitting and was also advocated in [|56[| , is 
to consider the ratio 

**<'■■*■> = ^T^T - {zL? - (43) 

Here, however, one needs a reliable signal for fairly large times t± + t 2 . Since the authors of 



36| accumulated 220 configurations they were able to apply this method successfully, which 
turned out to be consistent with the other ones. In view of our smaller statistical sample, 
we did not use the ratio R Z L(ti,t 2 ) to extract y^ tatlc . 

In order to determine the B parameter we computed the relevant 3-point correlator using 
the following expression 

K? s {t u t 2 ) = E (0|T{At 5 (f 1 ,-t 1 )aJ rtt (0)i4t 5 (f 2 ,t a )}|0) 

* i.*2> {Z 3 ) 2 E(ti+t 2 ) /do platt 

2M B ^ 1 

where i = L, R, S, N labels the four operators in the lattice effective theory, and A s (x, t) is 
the smeared axial current. In order to cancel the contributions from Z s and the exponen- 
tials in eq. (fP|) we consider suitable ratios of Kf s {ti ) t 2 ) and the 2-point correlator C SL (t), 
i.e. with the local operators always at the origin and the smearing performed at the same 
timeslices in both the numerator and denominator. 

Using the definition of Bb{h) in the continuum theory eq. and defining the ratio 



B°) , (44) 



Rf (ti,t 2 ) in the lattice effective theory as 

^ b * 2) = f^(|^ (45) 
then, provided t\, t 2 3> 0, the B parameter at scale is obtained from 

]T Z % (ZT tic r 2 Rf S (h, h) tl ^° B B (m b ), z = L,R,S,N (46) 



with the Zj's given in eq. 

In the computation of the ratio Rf s {t 1 ,t 2 ) on a periodic lattice we expect a signal for the 
correlator describing B° — B° mixing for t\ and t 2 on opposite halves of the lattice |57|], i.e. 



< tx < T/2, 
T/2 < t 2 < T, (47) 

where T = 48 is the length of our lattice in the time direction. In order to exploit time- 
reversal symmetry and thus to enhance the signal for the correlator, we compute the ratio 
Rf s (t 1 ,t 2 ) from 

i 2) ^ C SL i t l)+C SL {T _ tl) ^ C SL { t 2)+C SL {T _ t2)] ^ ) 

The correlators were calculated for timeslices 2 < t\ < 12 and 36 < t 2 < 46, which includes 
the entire region where one expects their asymptotic behaviour. 



3.3 Results for B 8 ^ 10 

In this subsection we present the results for the B parameter _Bg atlc using different smearing 
methods. For gauge-fixed configurations the four types of smearing defined in eqs. ( [T7| ) - (|20| ) 



were used, i.e. exponential (EXP), gaussian (GAU), cube (CUB) and double cube smearing 
(DCB). The results from the gauge- invariant Jacobi smearing algorithm are labelled (INV). 

We begin by describing the two methods we used to extract i?f^ atlc (/z) at renormalisation 
scale /i = m b from the ratios Rf s (ti,t 2 ),i = L, R, S, N. 

Method (a): The four ratios Rf s (t 1 ,t 2 ) are fitted individually to their asymptotic values 
Ri. The B parameter is then obtained through 

Bt tic (m b ) =J2 Z i (ZT ic )' 2 (49) 



with the factors Zi given in eq. (p7|). 



Method (b): Using the four ratios Rf s (ti, t 2 ), we define the B parameter S^ atlc (mfc; t\, t 2 ) 
for each set of timeslices (ti,t 2 ) 

5g atic (m 6 ; t u t 2 ) = Y,Z l (ZT ti T 2 R? S (ti,h), (50) 

i 

and fit Bg atlc (mb, t\, t 2 ) to a constant in suitably chosen intervals of (ti,t 2 )- 

The plateaux in the ratios Rf s (ti,t 2 ) can most conveniently be identified by fixing t\ at 
£1 — tf < T/2, and studying Rf s (ti,t 2 ) as a function of t 2 only. 

In order to illustrate method (a) we show in Figure |I] the plateaux for the four ratios Rf s (tf, t 2 ) 
for tf = 3, using Jacobi smearing at k\ = 0.14144. It is seen that a good signal is obtained 
for the four lattice operators, on the backward half of the lattice as expected. 

In Figure^] we show the signal for i?§ atlc (m&; tf, t 2 ) obtained using method (b) for both cube 
and Jacobi smearing. Despite the slightly shorter plateau for cube smearing which is also 
observed for all other smearing types in Coulomb gauge, the signal obtained in this fashion 
is also very clear. 

As one goes to smaller quark masses, the signals become noisier but are still of good quality, 
and the plateaux can easily be identified. For the other smearing types computed in the 
Coulomb gauge the picture is similar, and therefore we do not show the plots corresponding to 
Figs.[l| and|[ In general, B^ atlc (mb] tf, t 2 ) and the ratios Rf s show slightly larger statistical 
errors when Jacobi smearing is used, but, apart from that, the different smearing types give 
very similar results. 

We have also studied the behaviour of the plateaux for different values of tf. Using a larger 
value, e.g. tf — 4 leads to bigger errors in the ratios Rf s (tf,t 2 ) and Bg atlc (mb]tf,t 2 ), and 
the plateaux are shifted to smaller values in t 2 . On the other hand, using tf = 2 gives 
smaller statistical errors, but the plateaux are less flat which leads to higher x 2 /dof when 
fitting the ratios to a constant. We emphasise that the ratios are statistically compatible 
for tf = 2, 3, 4 for all smearing types considered. Hence, for our best estimates, using either 
method (a) or (b), we perform a simultaneous fit to the plateaux observed for tf = 3, 4. The 
results obtained by combining the plateaux for tf = 2, 3, 4 are quoted as a systematic error 
on our final value for B s ^ atlc . 

The values for B s ^ tlc {mb) extracted using methods (a) and (b) are entirely consistent. We 
have a slight preference for method (a): it gives better plateaux and offers more flexibility 
in the fitting procedure by ensuring that each of the four contributions are fitted in a region 
where the respective asymptotic behaviour has been reached. In the following therefore we 
base all our estimates on method (a). Correlations between different timeslices were taken 
into account in each fit. We did not attempt a simultaneous fit allowing for cross-correlations 
between different operators and timeslices, as systematic effects among the four operators 



2 



1 — 



R L (t f ,t 2 ) 



I I I 



X /dof = 0.38/4 



35 



40 



45 



R H (t f ,t 2 ) 



1 — 



_L — ^ 



X /dof = 0.70/4 



35 



40 



45 



0.0 



-0.5 — 



R 3 (t„t 2 ) 



-1.0 



-1.5 — 



X /dof = 0.41/2 



-2.0 



35 



40 



45 



R N (t t ,t 2 ) 



1 — 



— K light = 0.14144; smear = INV; t f = 3 
X 2 /dof = 0.49/3 



35 



40 



45 



Figure 1: The ratios R s L s {t f ,t 2 ) ) R s R s (t f ,t 2 ), R s s s (t f ,t 2 ) and Rff(t f ,t 2 ) for t f = 3 and 
Ki = 0.14144 using Jacobi smearing. The solid lines represent the fits over the respective 
time interval. 
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Figure 2: The quantity Bg atlc (mi,; tf, t^) defined in eq. (EDI) for cube smearing (left) and 
Jacobi smearing (right) at K\ = 0.14144 and tf = 3. 



could well be different. As was noted in ref. ||38|| , this could lead to an amplification of 
systematic errors in the fitted values. 

In Table |3T3| we show the results from fitting the ratios Rf s (ti, t%) to a constant for the three 
light hopping parameters and for all smearing types considered, and also list xVdof and the 
value for B s ^ atlc (mb) after renormalisation according to eq. fl49|) . A remarkable feature is the 



consistency of the results among all smearing types considered, which we take as evidence 
that the asymptotic behaviour has been reached. 

The values for the ratios Rl, Rr and Rn of the operators which mix due to the explicit chiral 
symmetry breaking induced by the Wilson term, are close to one, which is in accordance 
with the expectation from the vacuum insertion approximation. 

One notices only a weak dependence of the four ratios and of B s ^ atlc (mb) on the light quark 
mass. In fact, the results for Bg atlc (mb) are compatible with a completely flat chiral be- 
haviour within statistical errors as was already noted in |12[ . Assuming a linear dependence 
on the mass of the light quark, we can now extrapolate our results for B s ^ atlc (mb) to the chi- 



ral limit at K crit = 0.14315(2) or to the mass of the strange quark which, according to [15 |, 
is found at k s = 0.1419(1). Figure |] shows the chiral extrapolations of Bg atlc (m&) for cube 
and Jacobi smearing from both correlated and uncorrelated fits. Although the measured 
values appear to be almost perfectly linear as a function of the quark mass, the correlated 
extrapolation misses the points at smaller quark masses which might signal the use of a bad 
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Table 1: Results for the fits to the ratios of the four lattice operators using method (a). 
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Figure 3: Chiral extrapolations of _Bg atlc (m&) for cube smearing (left) and gauge-invariant 
smearing (right). The solid lines denote a correlated chiral extrapolation, whereas the uncor- 
rected fits are denoted by the dashed line. The extrapolated values from both procedures 
are shifted slightly in am q . 



fitting function. This results in a higher value for B^ atlc (m b ) in the chiral limit than that 
from the uncorrelated fit. The values for 5§ atlc (mb) from the two extrapolations agree well 
within errors, but except for Jacobi smearing the correlated fits have fairly large x 2 /dof. 

We quote our best estimates for _gg atlc (m b ) at /t crit and n s from the correlated chiral ex- 
trapolation. The measured values at the three light hopping parameters are shown with 



the extrapolated results in Table pT3| . In the following we will quote our best estimate from 
Jacobi smearing, which gave the cleanest chiral extrapolation. 

Combining the spread of values obtained from using different smearing functions and fitting 
intervals, extracting _gg atlc (m^) using method (b), and from performing uncorrelated fits into 
an estimate of systematic errors, we find 



B Bd (m b ) = 0.69 ^(stat) + _\ (syst) 



B Bs (m b ) = 0.71 ^(stat) ^(syst). 



(51) 
(52) 



At leading order, the renormalisation-group-invariant B parameter _Bg atlc is obtained from 



B s g Atlc {mb) according to eq. (|25f) for n/ = 5 active quark flavours via 



rpstatic 

B B 



1.476 B s ^ tic (m b ) 



(53) 
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Table 2: Correlated chiral extrapolations of B B atlc (mb) to K cr i t and k s for all smearing types 
considered. 



where we have taken A^ = 130 MeV in the expression for a s , eq. 



MS 



the relation between A^ and A^ given in |39| . We obtain 



B, 
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1.04 



+ 5 



(stat) ^ (syst) 



(stat) ^ (syst) 



in accordance with 

(54) 
(55) 



Within our errors we conclude that in the static approximation the matrix element of the 
four-fermi AB = 2 operator is indeed consistent with one. However, if the matching of matrix 
elements between full QCD and the lattice effective theory is performed using eqs. (|3l|) and 
(|34]), rather than eq. (|35|) , as was discussed in subsection the above values change to 



B, 



1.19 + _ 5 6 +3 2 and B B 



1.21 



5 * , respectively. Therefore we conclude that our 
and (|55|) are subject to a further 20 % uncertainty arising from 



best estimates in eqs. (| 
higher-order contributions to the renormalisation constants. A method demonstrating how 
to determine these factors non-perturbatively was discussed in [B3 



In previous lattice calculations of the B parameter using the static approximation [61 12 



the matching factors Zj, i = L, R, N, S were still unknown. The authors of obtained 
results for the ratio Rl, which are consistent with our findings (see TablepO). 



In two other simulations [|D| ^T|, propagating heavy Wilson quarks with masses around 
wicharm were used to compute the B parameter. In ref. JH]] results were quoted for the 
quantities B 1 ^ and B^, which correspond to our definitions of Rl, Rr- Extrapolating 
their results to the mass of the B meson, the authors of |4(| find B 1 ^ = 1.01 ± 0.06 ± 0.18, 

Rlatt 
B LR 



1.16 ± 0.01 ± 0.11, which again is in agreement with our results for Rl, Rr in 



Table Performing a similar extrapolation in the heavy quark mass, the authors of [^T| 



compute the renormalisation-group-invariant B parameter as Bs d = 1-16 ± 0.07, which is 
not incompatible with our result at infinite quark mass, given the additional perturbative 



uncertainty in Bs d and Bb s - We wish to stress that the calculation of the B parameter 
should be repeated with propagating heavy quarks using an 0(a)-improved action in order 
to study 1/tuq corrections to our result by analysing the mass dependence of B B . 



3.4 Results for f§ atic 

In this subsection we present our results for Jj^ tatlc extracted using different smearing func- 
tions. As outlined in subsection |3.2| , our best estimates are obtained by first fitting C s (t) 
to extract the binding energy E and (Z s ) 2 . The value of Z s is then combined with the ratio 
Z L /Z s obtained from a fit to C LS (t)/C ss (t), in order to extract Z L . 

It has been suggested that this method of determining Z L and subsequently f^ atlc potentially 



suffers from an incomplete isolation of the ground state |42|| . Failure to extract the ground 
state results in higher values for the binding energy and consequently in higher values for 
Z s and hence Z L . Therefore the authors of E2J followed a variational approach, based on 



smearing functions obtained using a relativistic quark model. The variational approach was 



also used by the authors of |43| who constructed the complete set of smearing functions 



allowed by the cubic group for a given lattice size. 



In a recent study by the APE Collaboration |I3| a number of checks for the isolation of the 



ground state without using the variational approach were presented. As will be shown later 
in this subsection, our results for y^ tatlc are entirely consistent with those in f44}| . 

Here, for the gauge- fixed configurations, in addition to exponential smearing (EXP), we also 
used a radially-excited exponential smearing function (EXP2S) defined by 

f(x, x') — \x — x'\ exp { — \x — x'\/r } (56) 

which is expected to have a considerable overlap with higher states. We then employed a 
variational approach to estimate the size of possible contamination of the correlators from 
the first excited state. Our values for E, Z s and Z L from 2-state fits were then compared 
to the results obtained from the other smearing functions using the procedure outlined in 



subsection 3.2. 



Following ref. we constructed a matrix correlator Cf. ■ (t) using the EXP and EXP2S 



smearing types as a 2 x 2 variational basis and determined the eigenvalues and -vectors of 
the generalised eigenvalue equation 

Cf(f + 1) vf = K(t+ 1, t) Cf (t) vf\ (57) 

For sufficiently large times t, the eigenvalues X a approach the eigenvalues of the transfer 
matrix, and therefore 

Xi(t + l,t) ^ e~ E (58) 



(59) 



where E, E* are the binding energies of the ground and first excited states respectively. 
Hence the quantity 

4ff(t) = log A^TM) (60) 



approaches the energy difference A£? = E* — E for sufficiently large i. In Table |Q| we show 
the values of 8 e g(t) as a function of i for all three values of Ki. It appears that 5 e ff(t) shows 
a plateau already for times around t = 2. Therefore we fix Ai? to be S e s{t = 2) and perform 
a constrained 3-parameter fit of the correlation function C (£), S = EXP according to 

C ss (t)^(Z s fe- Et {l + ^e-^} (61) 

where Z s * is the amplitude of the first excited state. It is this fitting form which we will from 
now on call a 2-state fit, whereas the usual 1-state fit is defined in eq. (|39|). It is possible in 
principle to use the eigenvector corresponding to the eigenvalue Ai to project the matrix 
correlator onto the approximate ground state. However, the resulting correlation function 
does not differ appreciably from the one computed using the usual IS exponential smearing 
function, and therefore we did not pursue this possibility further. 
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Table 3: The effective energy difference for the first three computed timeslices at all values 
of the hopping parameter of the light quark. 



In order to compare the results from the 1- and 2-state fits, we follow ref. |42| and plot the 
results for E and Z L obtained from 1-state fits to both EXP and EXP2S correlators over a 
time window t m - m , function of exp{— AE t m i n }. Keeping the length of the fitting 

window fixed at t max — t min = 6 and increasing t min allows one to extrapolate the results 
from the 1-state fits to t — oo. Therefore, as t min is increased, we expect that the results 
from 1-state fits converge to the value obtained from the 2-state fit performed over a large 
interval in t. 

In fact, as Figure^ shows, the results from the 1-state fit for the IS exponential smearing 
function (EXP) are in agreement with the 2-state fit already for £ min = 2. Furthermore, the 
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Figure 4: Values for the binding energy and Z L obtained from 2-state fits compared to 
1-state fits performed using a "sliding" window t min , t max , as a function of exp{— AE t m i n } 
for exponential smearing functions at Ki = 0.14144. 

results for the 1-state fits using the 2S exponential smearing function, which is supposed to 
have a poorer overlap onto the ground state, show indeed the expected stronger dependence 
on i min . We conclude that the fitting form eq. (|39|) applied to the IS exponentially smeared 
correlator is able to isolate the ground state correctly, provided the fitting interval is chosen 
suitably. Thus the 2-state fit merely serves to confirm the result obtained using the 1-state 
fit. This observation is further strengthened by comparing the fitted values for E, {Z s ) 2 and 
Z L for 1- and 2-state fits. The comparison is shown in Table 1X1 for one value of the light 
hopping parameter. 

The CUB and DCB smearing types are more problematic. Here, the plateaux for the C (t) 
correlators only start around t min = 6, 7 compared to t min = 4, 5 for EXP and INV smearing. 
The results from 1-state fits for the binding energy, however, are quite consistent with the 
results from EXP and INV, whereas the values for Z L are higher. Using AE determined 
from EXP smearing in a 2-state fit of the CUB and DCB correlators results in lower values 
for Z L but with a significant increase of the statistical errors. Single cube smearing (CUB) is 
particularly bad in this respect. One may suspect that the cube size was not tuned correctly 



in order to optimise the overlap of the operators. However, in [46] it was shown that sizes 
of ro = 4 and 6 gave substantially worse results than r = 5. 

At any rate, the values for Z L from all smearing methods differ by at most one to two stan- 
dard deviations, which is remarkably consistent, given the very different smearing functions 
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Table 4: The binding energy E, {Z s ) 2 , the ratio (Z S *) 2 /(Z S ) 2 between the first excited 
state and the ground state, and the final result for Z L for 1- and 2-state fits for exponential, 
and l-state fits for gauge-invariant smearing at Ki = 0.14144. Also shown are the respective 
fitting ranges. 



employed to enhance the signal of the ground state. 

We have thus established consistency between the results from l-state and 2-state exponen- 
tially smeared correlators, plus consistency among exponential and gauge-invariant smearing. 
We have also checked the stability of our results when directly fitting the correlator C (t) 
and computing Z L from J R(t) x Z L Z s . As was reported in ||46|| , results from exponential 
and gauge-invariant smearing are stable under the variation of the fitting procedure, whereas 
CUB and DCB smearing exhibit greater sensitivity to the method and fitting ranges em- 
ployed. 

The results for the binding energy E and Z L for all values of Ki are shown in Table |3.4| . Also 
shown are the extrapolated values at K crit and n s which were obtained assuming a linear 
dependence of E and Z L on the light quark mass. 

In the following we will take the results from the 2-state fits of the exponentially smeared 
correlators as our best estimate. Thereby we ensure that the more conservative choice of 
a larger statistical error encompasses all systematic variations in the final numbers from 
using gauge-invariant smearing and/or different fitting procedures. Thus, we do not quote 
an additional systematic error, and our final answer for Z L at n cr[t is 

Z L = 0.112 +g (62) 

7L 



which is in excellent agreement with ref. [36 in which Z = 0.111(6) is quoted. At the 
common value of K\ = 0.14144 in this work and ref. |[36|| , the values of E and Z L are consistent. 
Therefore we conclude that the small discrepancy between the binding energy obtained by 
APE and that in our earlier work [|l(J based on a subset of 20 configurations, has been 
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Table 5: Results for the binding energy E and Z L for all smearing types and values of K\. 
Also shown are the extrapolated values at K crit and k s . 



resolved. 

Using Z^ tatic = 0.78 and a^ 1 = 2.9(2) GeV we obtain 



yjtatk = 266 



L8 (stat) +J(syst)MeV 



-20 



(63) 



where the systematic error is due to the uncertainty in the lattice scale. Using the value of 
Z L at n s we obtain the ratio 

?jr- = 1.16 +3. (64) 

We can now compare our findings to other simulations. The direct comparison of f^ atlc [MeV] 
is however obscured by the different treatment of systematic effects such as the choice of 
^static an( j ^ e q Uan tity used to set the lattice scale. Therefore we choose to compare the 
results for Z L from simulations using the Wilson action [|12L |35|, EH, [43], |47], [42] and the 



0(a)-improved SW action [^, 44 1. Following a suggestion in ref. and assuming a scaling 
behaviour \ogZ L ~ log a and g~ 2 ~ log a, we plot \ogZ L as a function of f3 in Figure^. It 
is seen that the results (with the possible exception of ref. [[35|| ) are well aligned for {3 > 5.9, 
which supports the argument that scaling occurs in this region of j3. Furthermore there 
is consistency between the results coming from the variational approach (|43|, f|2| and this 
work) and those using the conventional approach |47], 36, |44 . 

The most striking observation is that, as far as Z L is concerned, there is practically no 
difference between the results obtained using the Wilson action and the SW action. A 
direct comparison was carried out by the APE Collaboration at j3 = 6.4 |Q and (3 = 6.0 



20| , and no difference within the statistical errors could be found. At first sight this may 



not seem surprising, since in the static theory improvement is performed in the light quark 
sector for which its effects on the mass spectrum were found to be small fl4|]. However, the 
renormalisation factor Z^ tlc at leading order in a s is quite different for the improved and 
unimproved action.^ In fact this difference amounts to an increase of 10-15% in the case of 
the SW action in the current range of (3 values. Consequently, collaborations working with 
the SW action quote relatively high values for f^ atlc in general, compared to those using the 
usual Wilson action. 

We conclude that at present the most severe systematic error in Jj^ tatlc is the uncertainty in the 
renormalisation factor Z^ tatlc . As in the case of the corresponding factors for B° — B° mixing, 
the perturbative estimate for this constant results in a large correction which signals that 
higher-order contributions may be important. A non-perturbative determination of 2'^ atlc 



using the method advocated in p3| for both the Wilson and the SW action is therefore of 
utmost importance. 



5 The one-loop expressions for Z s ^ atlc for the 0(a)-improved theory were computed independently by the 



authors of 
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Figure 5: Z L plotted logarithmically versus (3 using data from several simulations obtained 
using both the Wilson and the SW action. 



Systematic errors in J^ tatlc coming from uncertainties in the lattice scale will be further 
reduced once quantities that show good scaling behaviour, such as the 1S-1P splitting in 
charmonium or the hadronic scale Rq, become available for a wide range of (3. 



3.5 Phenomenological Implications 

We can now combine our best estimates from the previous two subsections into estimates 
for fsy/Bs- Using the results in eqs. fl51| ) and ( |63| ) we find 

fB d = 269 tf 5 (stat) +« (syst) MeV. (65) 

Using the values obtained after extrapolation to the strange quark mass we also quote the 
phenomenologically interesting ratios 

+ *(syst) (66) 
! 3 (syst) (67) 

!jM0, (68) 

where the systematic error is obtained from the spread of values using the systematic errors 
on B Bd , B Bs , as well as the result for Z L from a 1-state fit. 

The phenomenological implication of our results for f Bd , eq. (|63|), or f Bd \J B Bd , eq. (|65|), is, 
however, uncertain due to a number of systematic effects such as 

• the lack of an extrapolation to the continuum limit 

• large uncertainties in the values of the renormalisation constants Z^ atlc , Z L , . . . , Z N 

• the need to account for 0(AQ CD /m b ) corrections 

• quenching, i.e. neglecting the effects of quark loops. 



fB a a/ B Bs 
fs d \J B B d 

fk B B d 

fl B Ba M Bs 
fk B Bd M Bd 



1.16 +J(stat) 

1.34 +g(stat) 
1.37 ^g(stat) 



In ref. it was shown that the extrapolation of Jj^ tatlc to the continuum can yield a result 
below 200 MeV (albeit with a fairly large upper uncertainty). We have performed a tentative 
extrapolation, combining our result with the result of ref. [fHfl . The extrapolation gave a 
central value of f^ atlc ~ 250 MeV at zero lattice spacing for the SW action^. The difference 
between the two results is partly due to the fact that Z^ atlc is significantly smaller for the 
Wilson action than for the SW action as we mentioned before. 



3 The extrapolated value does not change if Rq is used to set the scale instead of the string tension. 



Lattice estimates for f B , especially in the static approximation, should therefore be treated 
with caution for phenomenological purposes. However, it is reasonable to assume that sys- 
tematic effects partly cancel in ratios such as f B J f Bd - in fact, as was shown in [ 4"2f , the 



a dependence of this ratio is compatible with zero. Therefore, in the following we illustrate 



the effect of our findings on the CKM matrix, using only the ratios in eqs. flBTf ) to (|68|), which 
are considered to be less afflicted with systematic effects. 

We focus on attempts to constrain the CKM parameters p and 77 in the standard Wolfenstein 
parametrisation. The B® — B® mixing parameter Xd is given by 

C 2 M 2 

^ = r B d f 2 B B Bd M Bd fj Bd y t f 2 (y t ) \V td V tb \ 2 (69) 

where r Bd is the B® lifetime, f\ Bd parametrises short- distance QCD corrections, and f 2 is a 
slowly varying function of y t = m^/M^. The current world average for Xd is [^] 

x d = 0.76 ± 0.06. (70) 

Mixing in the B° s — B® system can now be exploited in order to place constraints on the ratio 
\V ts \ 2 /\Vu\ 2 : 

x s r Bs f, Bs M Bs f 2 Bs B Bs \V ts \ 2 



x d r Bd f] Bd M Bd f 2 Bd B Bi \V td \ 2 ' 



(71) 



In this ratio the dependence on the top quark mass is cancelled, and one is left with an ex- 
pression involving only the CKM matrix elements plus SU(3)fl avour breaking terms. Assuming 
Vb s = VB d , and taking our estimate for f Bs B Bs M B J f Bd B Bd M Bd , we find 



— = (1.38 ± 0.17) (72) 
x d \Vtd\ 2 



where we have used r Bd = 1.53 ± 0.09 ps and t Bs = 1.54 ± 0.14 ps. This result is in good 
agreement with ref. [fy] where the proportionality factor is quoted as 1.25. 



Using the experimental result for Xd, we will now infer a value for the mixing parameter 
x s . This requires an estimate for the ratio |Vy 2 /|V^| 2 , which is usually obtained from 
global fits using the better-known CKM matrix elements as well as experimental data and 
other theoretical estimates as input. Various analyses of this kind have been presented in 

2 /IT/ 12 



50|, E% I52L pi 154], EB. In the standard Wolfenstein parametrisation the ratio \V t d\ l\V t 



f.s 



reads 



^ = A 2 (l-2p + p 2 + r ? 2 ) (73) 

I V ts\ 

where A = |V^ S | = 0.2205 ± 0.0018 | 56fl . The contraints on the CKM parameters p and 77 



depend crucially on the actual values of f Bd \jB Bd and Bk- 

In a recent study the authors have obtained values for p and r\ based on choosing 
Bk = 0.8 ± 0.2 (which is in agreement with recent lattice data |57|) and on the top quark 




Figure 6: The mixing parameter x s as a function of f Bd J B Bd for a fixed value of Bk = 0.8 
using our result eq. (^). The solid line follows the central values, whereas the dotted line 
represents the error band obtained from the errors on f Bs B Bg M B J f B B Bd M Bd and xj- 

mass of m t = 174 ± 16 GeV from CDF |[58|| . Using their values for p and rj and our estimate 
for f Bs B Bg M B J f Bd B Bd M Bd in eq. QB"%|), we plot in Figure the B® — B® mixing parameter x s 



as a function of f Bd \jB Bd . It is seen that values of f Bd yB Bd > 200 MeV result in practically 
unmeasurably large values of x s > 20. The current experimental lower bound is 

x s >9.0 (95%C.L.) (74) 



The error band in the figure is obtained from the errors on our value in eq. ( |68"D and on the 
experimental result for Xd- The errors on x s should, however, not be taken too seriously, 
since variations in Bk introduce large uncertainties in the ratio |Vt s | 2 /|Vt£z| 2 . 

We conclude this section by noting that the minimum x 2 i n the global fits in |52|| occurs at 



larger values of f Bd J B Bd as Bk is increased. This indicates that large values like f Bd yB Bd > 
200 MeV, as observed in lattice calculations in the static approximation, favour Bk — 1. 
However at present Bk is not available with good enough accuracy to provide further hints 



on the possible range of f Bd J B Bd . 



4 Mass Splittings 

In this section we report on our results for the B s — B^, — B, and B* — B mass splittings. 
The B* — B splitting receives particular attention since it first arises at order I/tuq and 
therefore serves to test the quality of the heavy quark mass expansion. For all the splittings 
computed we make a comparison with the results using propagating heavy quarks and with 
experimental data. 



4.1 The B s — Bd mass difference 

The B s — Bd mass splitting is obtained from the chiral behaviour of the binding energy E 
extracted from fits to the pseudoscalar 2-point function according to eq. fl3"9]) O]. Assuming 



a linear dependence of E on the light quark mass, we fit the chiral behaviour of E according 
to 

E ^= A+B Kl-it) (75) 

such that M Bg — M Bd is obtained from 

E{k s ) - E(K clit ) = M Bs - M Bd = B \(- - — ) . (76) 
The results for different smearing types, 1-state and 2-state fits, as well as correlated and 



uncorrelated chiral extrapolations, are shown in Table 4.1 





EXP 


CUB 


DCB 


INV 


M Bs - M Bd 


1-state 


2-state 


1-state 


1-state 


1-state 


cor 


0.029 +1 


0.030 + l 

— 5 


0.028 +4 A 

— 4 


0.029 +A A 

— 4 


0.029 +3 


unc 


0.027 + _\ 


0.027 


0.028 + l 

— 5 


0.030 + i 

— 


0.028 



Table 6: The mass difference M Bs — M Bd in lattice units for all smearing types, using both 
correlated and uncorrelated fits. For exponential smearing the results from 2-state fits are 
also shown. 



Taking the correlated value of the 2-state fit to the exponentially smeared correlator as our 
best estimate, and using a -1 = 2.9(2) GeV to convert into physical units we find 

M Ba - M Bd = 87 (stat) + J (syst) MeV (77) 



where the systematic error combines the spread of values obtained from the uncertainty in 
a -1 [GeV], using the 1-state result and performing an uncorrelated extrapolation. 



This result is in excellent agreement with ref. |4"2]| , where a value of 86 ± 12 _ g MeV is quoted 
as the continuum result. In a recent high-statistics simulation by the APE Collaboration 



44j using the SW action at (3 = 6.2, the M Ba — Mg j splitting was quoted as 58 ± 14 MeV. In 



general, APE's results for a range of (3 values |H[ seem somewhat lower than those reported 
This is partly due to the fact that the string tension was used in j44j to set the 



in 



lattice scale, giving a lower value (e.g. a 1 = 2.55 GeV at (3 = 6.2) than we use. Converting 
APE's result into lattice units, one finds aM Bs — aM Bd = 0.023(6) which is to be compared 
to our determination of aMg s — aM Bd = 0.030 _ , and thus the two simulations are not in 



disagreement. Table flT| contains a collection of results in physical units from simulations 
using the static approximation. 



M Bs - M Bd [MeV] 



^7 



+ 15 + 6 



-12 



58 ± 14 

70 ±10 
86 ± 12 

71 ± 13 
96 ±6 



-12 



+ 7 



-16 



Ref. 



this work 
PI 



in 



42 1 



El 



in 



Comments 



static, SW, (3 = 6.2 

static, SW, (3 = 6.2 
static, SW, (3 = 6.2 
static, Wilson, a = 

static, Wilson, (3 = 6.0 
experiment 



Table 7: Our value for the B s — Bd mass splitting in physical units compared to other 
simulations using the static approximation. Also shown is the experimental value. 



The result in eq. (|77j ) can now be compared to the results using propagating heavy quarks 



T0|] : extrapolating the pseudoscalar mass splitting Mp s — Mp d linearly in 1/Mp d either to 



M Pd = oo or to M Bd , M Dd one finds 



M Bs 


~M Bd 


= 84 


+ 14 
-12 


tgMeV, 


M Pd 


= oo 


(78) 


M Bs 


~M Bd 


= 93 


+ 12 
-12 


t^MeV, 


M Pd 


= M Bd 


(79) 


M Ds 


~ M Dd 


= 107 


+ 12 
-12 


!>eV, 


M Pd 


= M Dd . 


(80) 



The result at M Bd = oo is in excellent agreement with the static result in eq. (|77|). Fur- 
thermore, the value at Mp d = M Bd agrees very well with the experimental result of 96 ± 



6 MeV ||56| . The experimental value for Mp> g —Mp, d is 99.1±0.6 MeV |56|], which is compatible 
with our estimate. 

We conclude that for the M Bs —M Bd mass splitting it appears possible to interpolate between 



the static result and those obtained using propagating heavy quarks. From the behaviour 
of the splitting with 1/M P the size of Aq CD /mQ corrections is estimated at around 10 % at 
the mass of the B meson. 



4.2 The A b - B splitting 

In order to study the mass splitting of the A& and the B meson, we define a smeared 
interpolating field Af(x, t) according to 

Af (x , t) = e ijk bi(x, t) J2 f(x, a? ) (u\x\ t) C l5 d k (x', t)) (81) 

x' 

where /(af, x) is one of the smearing functions in eqs. fll3|) and ([17|)— (p0|) , and C is the charge 
conjugation matrix. In the above definition the spin of the baryon is carried by the heavy 
quark field b(x). 

We define correlators of the A& according to 

Cg(t) = E(Af(x,t)At S ( ))^i(^) 2 e-^ i (82) 
C L A s b (t) = ^(A L a (x,t)At S (0)) ^ Z^l*-^ (83) 

X 

where S = EXP, CUB, DCB, INV. We then obtain the Af, — B mass difference from an 
exponential fit to the following ratio of smeared-smeared (SS) correlators 

gfW E,(A»(g,t)At>)) m const e _ (EAi _ £)t 
C ss (t) £^f(2,t)4>)> 

where 

E Ab - E = M Ab - M Bd . (85) 

We used the same smearing functions in the numerator and denominator of the ratio in 
eq. (p4[) , although there is a priori no reason why one should do so. However, we found that 
the uncertainty in the ratio was dominated by the baryon correlator, and therefore we did 
not expect any gain in trying to optimise the signal using different smearing functions for 
the mesonic correlator. The ratios defined in eq. (|84]) gave short but clear plateaux in the 
range 9 < t < 11. 

The same procedure can of course be applied to the local-smeared (LS) correlators Cj^{t) 
and C LS (t). However, we observed that the effective mass plots for the ratio of LS correlators 
do not show clear plateaux. In addition, the fits of the correlators tend to give estimates for 
the splitting that are higher by up to two standard deviations, which further suggests that 
the ground state is not completely isolated in LS correlators. 



M Ab - M B 




EXP 


CUB 


DCB 


INV 


0.14144 
0.14226 
0.14262 


0.22 + _\ 
0.18 +3 
0.16 


0.23 + _\ 
0.19 
0.17 + l 

— 5 


0.22 + _\ 
0.18 +3 
0.14 + l 

— 


0.23 + _\ 
0.19 +2 
0.16 +3 


^crit 

X 2 /dof 


0.14 +3 
0.01 


0.17 +3 
0.18 


0.16 *\ 
0.78 


0.16 +3 
1.88 



Table 8: The — B mass splitting in lattice units at all three values of the light quark 
mass and extrapolated to the chiral limit. 



The ratio of correlators eq. ([84]) was fitted to a single exponential for 9 < t < 11 at all values 



of K\. In Table |4~2| we list our results in lattice units. 

Exponential smearing gave the cleanest signal at all values of K\. Assuming a linear depen- 
dence on the light quark mass, we extrapolated M Ab — Mb to the chiral limit. Again, the 
results from exponential smearing showed a very good linearity and consequently gave low 
X 2 /dof in the chiral fits (see Table |4~2"D . Furthermore, correlated and uncorrelated extrapo- 
lations gave almost the same central values. In contrast, the CUB, DCB and INV smearing 
types gave differing, though statistically consistent, results for correlated and uncorrelated 
fits. The x 2 /dof's are, however, larger and fairly high for gauge-invariant smearing. 

We therefore take our best estimate from the exponentially smeared correlators. In physical 
units we obtain 

M Ab - M Bd = 420 + 1 ™ (stat) ^(syst) MeV (86) 
with the systematic error coming from the uncertainty in a -1 [GeV]. Our value can be 



compared to other simulation results and the experimental number in Table p~2| . Comparing 
with the experimental value, it is seen that our new result is a marked improvement over a 



previous simulation in the static approximation [T3[ . 

In Figure [7| we plot our static result together with recent data obtained with propagating 



heavy quarks p9fl . Our value, compared to an earlier study |]13| , is in much better 
agreement with the mass behaviour of the results using propagating heavy quarks. In fact, 
an extrapolation in 1/Mp of these results to the static limit would be compatible with 
our value within the (relatively large) statistical errors. For further discussion of the mass 
behaviour, the reader is referred to 



M Ab - M Bd [MeV] 


Ref. 


Comments 


49n +100 +30 
_ 90 -30 


this work 


static, SW 


720 ±160 + 3 ° 




la 




static, Wilson 


QtJQ + 55 + 27 
60J -45 -26 




m 




prop., SW 


458 ± 144 ± 18 




m 




prop., Wilson 


362 ± 50 




& 




experiment 



Table 9: Our value for the A& — Bd mass splitting compared to other simulations using the 
Wilson action and/or propagating heavy quarks. Also shown is the experimental value. 
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Figure 7: The A& — B splitting in the static approximation (diamonds) compared to other 
simulations and the experimental values for the A^ — B and A c — D splittings. Only statistical 
errors are shown. 



4.3 The B* - B splitting 



Within the framework of large mass expansions, the B*—B mass splitting plays an important 
role, since it appears at order I/uiq; at lowest order (i.e. in the static approximation) 
Mb* = Mb- At order 1/tuq the splitting arises due to the spin-dependent, chromo magnetic 
correction term to the quark propagator 

Sl(x,0) = l -±-^5(x) £ fdrV^^a^F^^V^O) (87) 

i<j 

where Vg(t,r) and V$(t, 0) are defined according to eq. (10|), and Fij is a lattice definition 
of the field tensor. 



Following the discussion in ||T3| , |6l|| we compute the B* — B splitting from the insertion of 
Sp(x, 0) into the correlation function. The usual static correlator is given by 

C (t) = - £ (7475 S Q (x, 0) 7475 S t (0, x)) (88) 

x 

where Sq(x, 0) is defined in eq. (0), and Si(x,y) is the light quark propagator. In addition, 
we define 

C a {t) = -£ ( 7 475 ^(s, 0) 7475^(0, x)) (89) 
For large time separations, the ratio R a (t) = C a (t)/Co(t) shows a linear behaviour 

Kit) = %M A, + BJ. (90) 
The splitting M B , — M B is then given by the linear slope B a according to 

Ml*-Ml = zJ^B a . (91) 
where Z a is the renormalisation constant of the magnetic moment operator of the heavy 



quark |62|, |63|. As in the case of the renormalisation constant defined in eq. (|36|), we 
insert the reduced value of the quark self-energy into the expression given in |63|] . Using the 
"boosted" value of the gauge coupling in the numerical evaluation of Z a at one loop, we find 



63 



Z u = 1.76. (92) 

This is a very large correction, which suggests that higher-order contributions are likely to 
be important and highlights the necessity of a non-perturbative determination of Z a . 

In our simulation the SS correlator C a (t) was calculated using only gauge-invariant (INV) 
smearing. In the LS case, for which more smearing functions were used, the linear behaviour 
of R„(t) could not be established reliably. Thus, we cannot compare different smearing types 
for the B* — B splitting and therefore restrict the discussion to gauge-invariant smearing. 
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Figure 8: The ratio R a (t) for gauge- invariant smearing at Ki = 0.14144. The fit to eq. fl90D 
is shown as a solid line in the fitting interval 2 < t < 5 and continued as the dotted line for 
larger times. 

The ratio R a (t) was fitted to the functional form in eq. fl9"0|) for 2 < t < 5 at all three values 
of Ki. Figure [8] shows the signal at Ki = 0.14144 together with the fit. It appears that in 
addition to the linear behaviour of R a (t) for very small t, there is also a linear increase for 
7 < t < 11, albeit with much larger statistical errors. Fits to eq. ( |90"D in this time interval 
lead to values of B G which are larger by up to two standard deviations than those obtained 
using 2 < t < 5. The fits at larger times are, however, very sensitive to small variations in 
the fitting interval. We regard the two-sigma deviation at higher t as a correlated statistical 
effect, and believe that the asymptotic behaviour is already observed for small t. Of course, 
this must be confirmed in future simulations with higher statistics. 

The results for the values of the linear slope parameter B a from correlated fits are shown in 
Table f4~3[ together with the linearly extrapolated value at K CT a. Using uncorrelated fits gives 
essentially the same results. The values for B a increase slightly with decreasing light quark 



mass, as was already observed in [13], but are also consistent, within the statistical errors, 



with B a being independent of the light quark mass. 







0.14144 
0.14226 
0.14262 


0.0133 + _\ 
0.0137 +7 a 

— o 

0.0137 + l 

— o 


^crit 
K s 

X 2 /dof 


0.0143 +s a 

— o 

0.0136 + _\ 
0.54 



Table 10: The fitted linear slope B a of the ratio R a (t) in lattice units at all values of Ki, in 
the chiral limit and at the strange quark mass, extracted from the smeared-smeared (INV) 
correlator. 



Multiplying by Z a A/3 and converting into physical units we find 



M 2 B * -M 2 B = 0.281 ^(stat) +^(syst) GeV 2 (93) 



Ml* - Ml s = 0.268 ^3 (stat) ^(syst) GeV 2 (94) 
with the systematic error coming from the uncertainty in a -1 only. A comparison with 



experimental data and other simulations is made in Table |4~3| . Our result for the B* — B 
splitting is lower than the experimental value by almost a factor of 2. Also, contrary to the 
experimental observation, our estimate for the B* — B s splitting is lower than the one for 
B* — B. This can partly be accounted for by the opposite chiral behaviour of the splitting 
seen on the lattice. It is interesting to note, however, that both the experimental and 
lattice determinations of Ml* — Ml s yield a result that is compatible within errors with the 
corresponding value of Ml* — Ml. 

The use of the 0(a)-improved S W action for the light quark does not lead to a considerable 
increase in the splitting, as a comparison with the result of ref. |TB[ shows. This is not the 
case for propagating heavy quarks. When the SW action is used for both heavy and light 



quarks ||10|| , one obtains a value that is also about half of the experimental result. This 
is still an improvement over the case of propagating heavy Wilson quarks |13[] which gives 
a value about 10 times below the experimental result. As far as the hyperfine splitting is 
concerned, we therefore conclude that the main benefits of using the 0(a)-improved action 
are obtained in the case of relativistic heavy quarks. Given the large uncertainty in Z a , the 
Eichten expansion and propagating heavy quarks give comparable results if the SW action 
is employed. The discrepancy between the lattice and experimental results may, at least 
partially, be ascribed to quenching effects, as has been argued in ref. 



M|* - M 2 B [GeV 2 ] 


M% - Ml [GeV 2 ] 


Ref. 


Comments 


n 981 + 15 + 40 

- 281 -16 -37 

0.27 ±0.05 -0.07 
0.202 + 76 +29 

-84 -27 

0.488 ± 0.006 


026g +13 +38 

0.508 ± 0.028 


this work 




[!>6| 


static, SW 
static, Wilson 
prop., SW 
experiment 



Table 11: Our value for the B* — B and B* — B s mass splittings compared to other 
simulations and experiment. The value of Z a used in [^3| has been evaluated using the 
reduced value of the quark self-energy and the "boosted" gauge coupling. 



5 Summary and Conclusions 

In this paper we have reported on the results from an extensive study of weak matrix elements 
and the spectroscopy of heavy quark systems using the static approximation. A large part 
of our analysis was devoted to studying different types of smeared (extended) operators used 
in order to improve the signal/noise ratio and the isolation of the ground state. Although 
exponential or gauge-invariant smearing gave the best signal for most quantities, all the 
smearing functions gave remarkably consistent results. In addition, the variational approach 
employed in the determination of f^ atlc demonstrated the compatibility of results obtained 
using this more refined fitting procedure with those from the usual single-exponential fits. 
Thus we are confident that we correctly isolate matrix elements and spectroscopy data from 
the ground state contribution of suitable correlators. 

We obtain a good signal for the various four-fermi operators relevant for B° — B° mix- 
ing. Our estimate for Bb in the static approximation is in agreement with its value in 
the vacuum insertion approximation. Regarding we note that our determination of 
Z L = /^ tatlc ^2/ M B I Z s 2 atlc is consistent with other simulations. 

Among the systematic errors present in this simulation, the most important (apart from 
quenching) are due to uncertainties in the renormalisation constants relating the matrix 
elements on the lattice to their continuum counterparts. These systematic effects manifest 
themselves most severely in our estimate for the B parameter, and in the case of f B ta,tlc , 
where there is practically no difference in Z L for the Wilson and the SW actions, yet the 
corresponding values of Z s ^ tlc differ by about 10-15 %. Also, the large value of Z a in eq. fl92"|) 
implies that higher-order contributions may be important in the perturbative evaluation of 
this constant. 

Our results for the B s — Bd and — B splittings compare very well with experimental 
estimates, although the statistical errors, especially for the A& — B splitting, are still large. 



The B* — B splitting obtained from a I/tuq correction to the static limit, however, does 
not agree with experiment. Using the SW action for the light quarks does not lead to a 
significant increase in the lattice estimate of M B , — M B . Future simulations using dynamical 
quarks may reveal whether the discrepancy between the lattice and experimental results is 
due to quenching effects. 

The static approximation, in conjunction with a refined numerical analysis, remains a valu- 
able tool in lattice studies of heavy quark systems. It plays the crucial role of guiding the 
extrapolation of results obtained using propagating heavy quarks to the mass of the b quark, 
by providing direct information at infinite quark mass. 

In future, one should concentrate on the analysis of systematic effects such as non-perturbative 
determinations of the renormalisation constants. In the case of Bb it would be highly de- 
sirable to repeat the calculation for propagating heavy quarks, preferably with an improved 
fermion action, in order to study the mass dependence. 

Note added: after completion of this work, we received several papers, [p5|H|69|l, presenting 
results for the B parameter |65|, |66|, ^7|, fs |65|, |69[| , the mass splittings Mb s — MB d , 

M B — M B * [0, B K |5B| and discussing phenomenological implications |57| . The conclusions 
of this paper remain unaltered, though, since the reported numbers are in agreement with 
our findings. 
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